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^ ' Abstract 
O , 

' We discuss the transition from a metal to charge or spin insulating phases 



characterized by the opening of a gap in the charge or spin excitation spec- 
tra, respectively. These transitions are addressed within the context of two 



^ ' exactly solvable Hubbard and tJ chains with long range, 1/r hopping. We 

o ■ 

^ ' discuss the specific heat, compressibility, and magnetic susceptibility of these 

^ \ models as a function of temperature, band filling, and interaction strength. 

I We then use conformal field theory techniques to extract ground state corre- 

lation functions. Finally, by employing the g-ology analysis we show that the 
charge insulator transition is accompanied by an infinite discontinuity in the 
Drude weight of the electrical conductivity. While the magnetic properties 
of these models reflect the genuine features of strongly correlated electron 
systems, the charge transport properties, especially near the Mott-Hubbard 

transition, display a non-generic behavior. 
PACS1993: 71.27.-ha, 71.30.-hh, 05.30.Fk 
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I. INTRODUCTION 

RecentlyS we introduced a Hubbard-like model which describes spin-1/2 Fermions on a 
chain of L sites, hopping with long-range amplitude, ti^rn = it{—iy~"^[d{l — 'm)]^^. Here, 
d{l — m) = (L/tt) sin[7r(/ — m)/L] is the chord distance between sites / and m on the chain 
closed into a ring (the lattice spacing a is set to unity). In the thermodynamic limit, 
L — s> oo, and for a fixed distance, (/ — m), the purely imaginary hopping ti^m = t*m i becomes 
ii,m —^it/{l — m), leading to the Hubbard Hamiltonian with "1/r-hopping" , 

H = f + UD= hmctc^^^ + uj^hi^^hi,^ . (1) 

l^m=l,a 1=1 

Here U is the strength of the usual local Hubbardi interaction. For even L we choose 
antiperiodic boundary conditions, so that the dispersion relation is linear in wave vector, 
namely, e{k) = tk for k = A{m + 1/2) (A = 27r/L, m = -L/2, ... ,L/2 - 1). For f/ = 
the Fermi sea is the ground state with all /c-states from k = — tt to kp = TT{n — 1 ) fille J, 
where n = {N^ + N^)/L is the total particle density (the "filling"). Note that the kinetic 
and potential energy operators are respectively odd and even under partity, and thus parity 
is, in general, not a good quantum number in our model. 

In the large U limit the Hamiltonian (|^) becomes the 1/r-tJ model including "pair- 
hopping" terms {Pd=o projects onto the subspace with no double occupancies)^: 



^tJ = pD=o[f + E 



- E ^EMqtctA-.'C_'}^D.o . (2) 

As usual, hi = hi^^ + n;^|, and is the spin-1/2 vector operator {S^' = c^|C; |, Sf = c'^'^iCi^^ 
Sf = {hi^i — hi^i)/2). At half filling, where the first and third terms in (Q) vanish, the 
model (^ reduces to the (parity symmetric) (l/r)^-Heisenberg model introduced by Haldane 
and Shastryi, with J = At^/U. 

In Ref. |1] we conjectured the full excitation spectrum and associated degeneracies, and 
were thus able to calculate the free energy. The exact solution allowed us to identify two zero 



temperature (T = 0) phase transitions: the first was a Mott-Hubbardi'i metal-to-charge- 
insulator transition (MCIT) in the half-filled 1/r-Hubbard model. This was signaled by the 
opening of a charge gap for f/ > Uc = W where W = 27it is the electron bandwidth. The 
second, a metal-to-spin insulator transition (MSIT) was associated with the opening of a 
spin gap in the 1/r-tJ model (@) for J > Jc = 2W/ [(1 - n)n^]. 

In this paper we clarify the nature of these transitions and attempt to separate the fea- 
tures due to the special form of the dispersion in (|I|) from more generic properties which 
may be expected to survive for other dispersions and/or in higher dimensions. The plan 
of the paper is as follows: we start in section || by recalling the form of our exact solu- 
tion for the spectrum and ground state energies, from which we identify the location of 
the MCIT and MSIT. In section |T| we discuss the thermodynamic properties of our two 



models, namely the specific heat, compressibility, and magnetic susceptibility. In sec. |V 
we use conformal field theory techniques to extract the long-range behavior of ground state 
correlation functions. Sec. contains a discussion of the connection with the "g-ology" 
approach to one-dimensional systems, which we exemplify by calculating the Drude weight 
in the frequency dependent electrical conductivity for the 1/r-Hubbard model. Finally, a 
summary and conclusions are presented in sec. |VI| . 

II. GROUND STATE PROPERTIES AND EXCITATION SPECTRA 

We recall the effective Hamiltonian, introduced in ref. |l|, 

-TT<K<TT I O" ) 

which describes spin (s^co-; = ±1/2, Cjc = 0) and charge {dic,ejc; Sjc = 0, C^^ = ±1/2) 
degrees of freedom in an occupation number representation with a hard core constraint, 
^an-^ca + ^K + ^l = 1 for each /C. Furthermore, /i^,, = (t/C)/2-/i^, /i^ = -{t}C)/2-2fi + U, 



h"^ = -(t/C) /2, and J/c = [t{2IC - A) - U + ^{27Tt)'^ ± f/^ - 2tf/(2/C - A) J /2 > 0. The 
chemical potential in the presence of an external magnetic field is given by = /i — 
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a{gfiB'Ho)/2. In the following we will set fiB = 1, and g = 2. To be precise, we also restrict 
ourselves to t > and U > — 27rt in which case we were allowed to identify /Cmin — A = /Cmax 
because Jjc vanishes for JC = /Cmin = — A(L — l)/2. 

Note that, formally, the entire spectrum appears to display spin-charge separation at all 
energies in the sense that the Hamiltonian splits up into independent spin and charge contri- 
butions ("strong" spin-charge separation). In reality, however, spin and charge excitations 
are coupled by the constraint, J2a ^a: o- + + ^k; — 1? spin-charge separation only oc- 
curs at sufficiently low energies/temperatures, where spin and charge excitations contribute 
independently to various physical properties. 

As already discussed above, we will concentrate on the physics of two special limits: 
(i) the l/r-Hubbard model in the vicinity of half filling, < 1, and (ii) the 1/r-tJ model 
for n < 1. The latter is obtained by taking the limit U —>■ oo oi which projects out all 
double occupancies (/i^ —>■ oo). With the help of the completeness constraint one arrives at 
the tJ effective Hamiltonian, 

H^! = Hf + Hf= Y: [(-t/CK - J^^^_a Ai] ' (4) 

— 7r</C<7r 

where, for J = At^/U ^ 1, the exchange coupling, J/c^^U/t), reduces to J/c = 
{J/4:) TT^ — (/C — . Below, we will adopt the usual standpoint for the tJ modelB, 

and treat J as an independent parameter. 

A. 1/r-Hubbard Model 

From eq. we can immediately extract the form of eigenstates in an occupation number 
representation in terms of the effective spin {sic,a) and charge (e^c, d/c) degrees of freedom. 
For example, the ground state is expressed symbolically as 



ground state: [Ti] • • • [Ti] 



o . . . o 

/CF=T(2n-l) 



where t, |, •, and o represent Sjc^^, sic^i, die, e^:, respectively. We brace those pairs at /C — A 
and JC which contribute an interaction Jjc i. The ground state can be regarded as a short- 



range RVB state in /C-space obtained by filling K, states with [flj-pairs from K, = — 7r + A/2 
to /C = ICf — A/2, where JCp = vr(2n — 1) . Although the real-space structure of the 
ground state wave function is not immediately apparent, it should be clear, however, that 
it also consists of long-range (overlapping) singlet pairs. In fact, the recently constructed 
wave functions for the U = oo limit (the "l/r-t-model")i are genuine RVB stateJi^l of the 
Gutzwiller-Jastrow type. 

The corresponding ground state energy density, eo{n), can easily be obtained by apply- 
ing the effective Hamiltonian (^) to the symbolic ground state wave function given above: 
eo(n) = (l/27r) j'^^ dJC [J(/C)/2] + (l/27r) /^^ rf/C(-t/C), where the factor 1/2 in the first term 
takes into account the fact that only every second /C value contributes to the first integral. 
A simple calculation then gives 

Un-W{l-n)n 1 



eoin < 1 



{W + Uf - i{W + Uf - mun^''^ 



(5) 



4 2AWU L 

where W = 2Tit is the bandwidth. Using particle- hole symmetry0 we obtain coin > 1) = 
eo(2 — n) + U{n — 1), and correspondingly, for the chemical potential at zero temperature, 
/i(n < 1) = deoin < l)/dn, /i(n > 1) = f/ — yu(2 — n). From (^ one can easily see that 
A/ic = fJ'+{n = 1) — /i„(n = 1) = (lim„^i+ — limji^i-)/i(n) becomes finite for U > Uc = W: 
= U — Uc- Equivalently, the charge compressibility, k, = dn/dfi at n = 1 vanishes for U 
above Uc, as expected at a Mott metal-to-charge- insulator transition (MCIT). 

We are now in the position to identify the low-lying spin and charge excitations which 
we discuss separately for the half-filled and less than half-filled cases: 

1. Half- filling: n = 1 

We begin with the spin excitations. The four degenerate 5* = 0, S* = 1 lowest-lying spin 
excitations are represented by 



spin-exc: [Ti] • • • [Ti] o" 



[n]...[n] ^' 



[Ti] • • • [Ti] 

fC2 
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For /C2 = A^i + A there is a triplet excitation only. In the thermodynamic limit the excitation 
energy is given by 6Es{ICi, IC2) = [^(/Ci) + -^(^^2)] /2, and we may thus identify the lowest- 
lying spin excitation with two spinous. They are characterized as spin-1/2 objects which 
always come in pairs, and are separated in /C-space by an even multiple of A. We thus 
rescale /C = 2/C' to retain the proper spacing A = 27r/L of /C' values. The spinous then have 
the dispersion relation, 

e,(/C) = J(2/C)/2 = (^^W^ + U^- mUK/'K + 2W1C/'k - t/^ /4 > ; |/C|<7r/2, (6) 

which we depict in fig. 0. The spin excitations are always gapless at /C = — 7r/2, and at 
/C = 7r/2 for U /W > 1. Their corresponding velocities at /C = ±7r/2 can now be calculated 
from fs(/C) = (9es(/C)/(9/C as 

vf^Vs{JC = -n/2) = jjj^—^ for all U/W (7a) 
v^^v.,ilC = n/2) = -jjj^-^ for U/W > 1 . (7b) 

Here, vp = t is the Fermi velocity of the bare particles. The density of states 
for spin excitations is then calculated as Ds{E) = {1/ L)Y,--k/2<k<-k/2^ i.^ ~ ^s{]^)) = 
l/(27r) Xr{/2 S {E — es(/C)). With the help of eqs. (|^), for low energies we obtain 



DsiE^O) 



1/ (27rt;f ) for U/W < 1 

1/ (27rt;f ) + 1/ [27t\v^\) for [//W^ > 1 



We next turn to the charge excitations. We restrict ourselves to a fixed particle number 
(C^ = 0) and we only obtain two of the four C = 0, C = 1 lowest-lying states. These charge 
excitations are represented by 

charge-exc: [Ti] . . . [Ti] o [Ti] • • • [Ti] • [Ti] • • • [Ti] ■ 

For /C2 = /Ci + A there is a low lying triplet excitation [C = 1, = 0), with the singlet 
excitation (C = 0, = 0) always at high energy (anti-bound state). Again, we may identify 
the lowest-lying charge excitation with two chargeons. In the half-filled case we always get 
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a pair of a holon and a douhlon which are separated by an even multiple of A = 27r/L. 
Rescaling /C = 2/C' leads to the chargeon dispersion relation, 



e^(/C) = L/W^2 + f/2 - 4W/C/7r - 2H^/C/7r + /4 > ; |/C| < 7r/2 , (9) 



as depicted in fig. 0. The charge excitations are gapless at /C = 7r/2 for [//VT < 1, with a 
velocity given by 

^c=- ^^^l^ for U/W<1 . (10) 

For U > W a. gap, 2ec(7r/2) = U—W = A/ic, opens in the charge spectrum. Correspondingly, 
the density of states for charge excitations at low energies takes the form 

f 1/ f27r|t;,^|) for U/W < 1 
D,iE^O) = \ ^ ^ . (11) 

[ O [exp {-E/Ai2^)] for U/W > 1 

At this point a peculiarity of our model becomes evident: the charge velocity increases as 
a function of U /W and eventually diverges ai U = Uc = W. The effective charge mass m* 
which is connected to the density of states or the velocities by m*/m = Dc{0)/ Dq{0) = 
vf/\v^\ decreases as a function of U/W. This is in contrast to the Brinkman-Rice scenario 
for the Mott-Hubbard transition in that case m*^/m increases and charge excitations tend 
to localize close to the transitionllllli. Indeed, for the Hubbard model with nearest-neighbor 
hopping (cosine dispersion) m*{n,U) diverges at the MCIT for U > when the transition 
is approached from below half- filling (n ^ V 



2. Less than half-filling: n < 1 

In this case there are only right-moving spinous with velocity v^, eq. (|7aD. The spinous 
are further restricted to —it/ 2 < K, < Kip /I with ICp = n {2n — V) . The lowest-lying charge 
excitations are now given by holons alone while the doublons are always gaped: 



charge-exc: [Ti] • • • [Ti] ° 



K 



[Ti] • • • [T 



i]o... 
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Although holons need not come in pairs, they are still separated by an even multiple of A 
in /C-space. Rescaling again we find for the charge excitation 



e^(/C) = e,(/C)-e,(/C^/2) + t(/CF-2/C) - ti /2 < IC < Kp 12 . (12) 



The holons are gapless at /C = /Cf/2, and the corresponding velocity, = deuiJC) / dlC 
leads to 



= -t;^ I 1 + ^ = 1 . (13) 

'{W + Uf - AWUn 



In turn, this implies that the density of states for charge excitations at low energies still 



remains 



3. Ground-state compressibility and magnetic susceptibility 

The results for the ground state energy density eo(^), eq. (||), allow us to calculate the 
chemical potential /x = deo{n)/dn and the T = compressibility k = dn/dfi. By turning on 
a small external magnetic field Hq we can also obtain eo{m,n) (see Appendix eq. ([A^)), 
and the magnetic susceptibility, x = dm/dHo, where the magnetization density, m, is related 
to Ho by Ho = deo{m,n)/dm. We summarize our results as follows: 

(i) n < 1 or {n = 1, U/W < 1): 



U - W{1 - 2n) - J{W + Uy - AWUn 
/i = (15a) 

1 2 
^ = = (15b) 

-^Ivcl W{l + U/^{W + Uy - AWUn) 
m = xHo (15c) 
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(ii) n = 1, U/W > 1: 

vv 

fx{n = 1") = — (16a) 
W 

^(n=l+) = U (16b) 

W (16c) 

(16d) 
(16e) 

7t\v^\ 

The system is an incompressible charge insulator at half-filling and U > W. Note that in 





w 


1") = 






Y 


1^) = 


u - 


A/ic = 


u - 


K = 





m = 


X^o 




1 


X = 





the limit n — >■ 1~ the compressibility stays finite for U > W (for U = W, k = 4-\/l — n/W 
for n ^ 1). This is, nevertheless, consistent with eq. ( |16d|) since, in contrast with the 
situation of the usual Hubbard mo in our case the function n(/i) is not differentiable 

at fi{n = 1^) and fi{n = 1+). 

We further note that eqs. ( |15b| ), ( |15cl| ), and (|TB|) are the generic behavior expected 



of a "Luttinger-Liquid"ll3, a point which we exploit in our discussion of the ground-state 
correlation functions (sec. [TVI ) and the Drude conductivity (sec. |V]). 

B. 1/r-tJ Model 

Since the case of half-filling has already been considered by Haldan(S, here we will be 
mainly interested in the less that half-filled situation. From the effective t J Hamiltonian (^) 
it is easy to see that the hole kinetic energy favors all particles to be as close to /C = — tt 
as possible while the exchange interaction J^q = [J / A) tt^ — (/C — A/2)^ tries to distribute 
the particles symmetrically around /C = 0. We thus expect a transition at some critical Jc 
in this case. 

In general, the ground state can be represented as 



ground state: o . . . 



[Ti] • • • [Ti] 



o . . . o 

Ki2 
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where /C2 = A^i + 2?™ has to be determined from the minimization of the ground state 
energy. We find (see appendix ^) 

Wn{l-n) J7r2n2(3-2n) 



eo(n, J < Jc 
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12 

y - 1] n{3 -2n) + 3-3n + n' 

Jc / 



with /Ci = — TT and 

eo(n, J > Jc 



j7T^n{3 - n^) 



24 



24J 



(3-n2)(^-l)+2(3-3n + 



J2 



72 



(17a) 



(17b) 



with /Ci = -77 (n + 2W^/(7r^J)) > -vr. Here, = 2iy/(7r^(l - n)) is the critical couphng. 
Note that Jc is proportional to 1/(1 — n), i.e., it decreases with increasing (hole) doping. 
Below we discuss the two cases, J < Jc and J > Jc, separately: 



1. J <Jr 



For J < Jc we can make use of the 1/r-Hubbard model results (sec. [II A|) . In particular, 
es(/C) = J((7r/2)^ — K?)/2 for — 7r/2 < /C < /Cir/2, and the spinon velocity at /C = — 7r/2 is 



j7r/2 . 



It is amusing to note that this is precisely the result for the Heisenberg-chain with nearest- 
neighbor interaction^. This is consistent with the fact that the Gutzwiller projected Fermi 
sea \ipo) = -PD=o|Fermi-sea)@, the ground state wave function of the 1/r-tJ model, is also 
an excellent trial state for the Heisenberg chain with nearest neighbor interactionSH, as 
well as for the nearest neighbor supersymmetric tJ model (J = 2t)ii. As in sec. |II A| the 
corresponding low energy density of states is given by 

D,iE^O) = ^ for J<Jc . (19) 

For < J < Jc, the holon velocity is calculated as 
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^- {2t + Jn{2n - l)/2) = -(7r/2) [(1 - n) {Jc - J) + nJ] . (20) 

Note that the hmit J ^ is pecuhar: precisely at J — we have a free gas of holons with 
e/i(/C, J — 0) — —2t{K, — K,f/2) for — 7r/2 < K, < /Cir/2, but with allowed /C-values spaced 
by A/2. On the other hand, the limit J — » only gives half oi the excitations since the 
/C values are now spaced by A rather than A/2. This is because, for J > 0, half of the 
J = excitations develop a gap, corresponding to the energy required to break a spinon 
pair, J(/Ci?) > 0. Consequently, the low energy density of states for charge excitations, 

L>e(£;^0) = ^^ for 0<J<J, , (21) 

is only half as big for J —> as for J = 0. 



2. J > Jc 



For J > Jc the lowest-lying spin excitation can be represented as 

[Ti] ■ ■ ■ [Ti] o o 



spin-exc: o . 


. o a 


o . 


.oa' 


o . 


. o 






K. 




JC' 





ICi 



o . . . o 



IC2 



In the present case nothing can be said about the ground state wave function in terms of the 
original Fermions, since the link to the 1/r-Hubbard model can no longer be made. Again, 
the spinous come in pairs but they can now have arbitrary separation in the region outside of 
/Ci < /C, /C' < /C2. A single spinon has the excitation energy es(/C) = t{}C — IC2) + J(/C2)/2 = 
t{lC — /Ci) + J(/Ci)/2 and thus, a gap proportional to J — Jc opens in spinon spectrum: 



A//5 = 2e^(-7r) 
7r2(l-n) 



(1-7) [J(l + n)-Je(l-n)] 



(22) 



J — Jc then corresponds to the onset of a metal-to-spin-insulator transition (MSIT). 
For J > Jc a charge excitation is represented by 



charge-exc: o . . . o [t i] . . . [U] o [Ti] ■ - ■ [U] 



o . . . o 



IC2 
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and, after rescaling /C, the excitation energy for /Ci/2 < K, < /C2/2 is eh,(/C) = t{K,2 — 2/C) + 
(J(2/C) — J(/C2)) /2. The corresponding velocities are 

= -v^ = j7m/2 (23) 

at /C = /Ci/2 and /C = /C2/2, respectively. We see that, for low energies, the gapless exci- 
tations are parity- symmetric although the underlying Hamiltonian itself does not conserve 
parity. The low energy density of states for charge excitations, 

D,iE^O) = -^ + -^ = -^ for J>J, . (24) 
Zirvt^ ZT^wtf] Jn^n 



doubles at the MSIT. 



3. Ground-state compressibility and magnetic susceptibility 



The T = chemical potential, compressibility, magnetization, and magnetic susceptibi- 
lity can be calculated from eqs. (p!7|), (|A5|) , ( |A8|) , and ( |A9|) of appendix ^ . As in the case 
of the 1/r-Hubbard model we summarize our results (Jc = 2W/ (7r^(l — n))): 

(i) J < Jc- 

W{1 - 2n) Jir^nil - n) 



(ii) J> Jc 



TT 



'1 - n)Jr^ 



4 



2n ( ^ - 1 



n\vj^\ TT^ [{1 - n) {J^ - J) + nJ] 
m = xHq 

1 2 



X 



2n^J 8 
SJ 



(l-n^) 



J2 



(25a) 

(25b) 
(25c) 
(25d) 



(26a) 



12 



K 



1 1 

+ ^ 



4 



TCV 



R 



Ml 



n)] 



m 



for Hq < ni 

X = . 



(26b) 

(26c) 
(26d) 

(26e) 



The system is a spin insulator for J/Jc > 1- Note that the compressibihty k, doubles at 
the transition because another low-energy charge excitation replaces the spin mode which 
becomes gaped at J = Jc. 

III. THERMODYNAMICAL PROPERTIES 

As noted in Ref. |l], the effective Hamiltonian, (^, is equivalent to a classical Ashkin- 
Teller type model in the presence of an inhomogeneous "magnetic field"ii, with periodic 
boundary conditions. The thermodynamic properties of our system can then be calculated 
by transfer-matrix technique^l. We define the abbreviations Sjc^a = exp (^—Ph'l^^^^, Djc = 
exp (^—(3h'^^, Eic = exp (— /3/i^), and Pjc = exp {—(3Jjc) where (3 = l/fc^T {ks = 1) is the 
inverse temperature. The transfer matrix between sites /C — A and /C, 







Dk 


Ek 


'S'/c,T 




Dk 


Ek 


Sic,-] 




Dk 


EkPk 






Dk 


Ek 



I 



has two vanishing eigenvalues; the remaining two eigenvalues are given by 



Xk±x\X^ 



Sk^Sk^i{\-Pk')^DkEk{\-Pk) 



4 



where we introduced the abbreviation Xk = Sk,i + Sk,i + Dk + Ek- The partition function 
for all chemical potentials /x, magnetic fields TCq, interaction strengths U/W, and tempera- 
tures T, 
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^L = [ n + n > (27) 

leads, in the thermodynamic hmit {N^,L — > oo, = N„/L = fixed), to the free energy 
density f{fi,no,U/W,T) = Iurl^^ [- lnZi/(/5L)], 

1 die 

f{^,no,U/W,T) = -- -^lnA^+ . (28) 

P J-n ZTl 



Here we made use of the fact that Ax; < Xjc^ This general form (28) remains valid 
for the 1/r-tJ model (^ if we set Die = 0. From eq. (^ it is easy to see that at half- 
filling the problem is equivalent to an Ising model on a ring with nearest-neighbor coupling, 
J/c = (J/4)(7r^ — (/C — A/2)^), and one thus recovers all of Haldane's results for the spin- 
chain0. 

A. 1/r-Hubbard model at half filling 

First note that at half-filling (n = 1, i.e.0, yu(T) = U/2) without external field the 
spectrum is completely specified in terms of independent spin and charge excitations. In 
fact, the free-energy in that case has the simple spin-charge separated form 

f{n = 1, U/W, T) = -U/2 + Coin = 1) 

die {In [1 + exp (-/?e,(/C))] + In [1 + exp (-/?e,(/C))]} (29) 

J-n/2 



27Tp J-n/2 

where the dispersion relations for the spin ("up" and "down" spinous) and charge (holons 
and doublons) excitations were already given in eqs. (^, (|^) (see figure |^). Note that the 
rescaling of /C-values gives an additional factor of two in front of the integral such that ( [29| ) 
corresponds to the result of two independent free fermion systems for charge and spin. 

1. Specific heat 

The spinon and chargeon densities of states, Ds,c{E) = l/(27r) J'^^J'^^ dlC6 {E — es^c(^)), 
are shown in figure ^ and 0b, respectively. Already before the MCIT the density of states 
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for the spinons develops a van-Hove singularity at f/ = W/2. This could be interpretated as 
the formation of the lower (spinon) and upper (chargeon) Hubbard band which are, however, 
not yet separated. 

Given the densities of states, the internal energy density u{n = 1, T) can be obtained from 
uin = 1, T) = eo(n = 1) + 2 /o~ dE [DsiE) + D,{E)] Ef{E), where /(E) = [exp{(3E) + 1]-' 
is the Fermi-Dirac distribution function. The specific heat = —P'^du/dp is then given by 

cJn = 1, T) = 4T / dx (DJ2xT) + DJ2xT)) — — . (30) 
JO V cosh X J 

At low temperatures this reduces to c^(n = 1, T — > 0) = 7(72 = 1)T with 

7(n = l) = y[D,(0) + D,(0)] . (31) 



This Luttinger Liquid relation!!^, which remains valid for all fillings, will allow us to identify 
the conformal charge of the conformal field theory that governs the low energy behavior of 



our model (see sec. ^V]) . 

The behavior of the Sommerfeld coefficient 7 as a function of the interaction parameter, 
U/W, can be extracted from eqs. (0), (|), ([ig), and (|ll]): 



1 for U/W < 1 
7(n = 1) = -Do(0) <; (32) 

' U/W for U/W > 1 



where -Do(O) = 2/W is the non-interacting density of states. Note that 'y{n = 1) is 
unrenormalized below the MCIT. The reason is that the increase in the density of states 
for spin excitations is exactly compensated by the decrease in the density of states for the 
charge excitations. This precise cancellation of the two contributions is a peculiarity of our 
model and does not occur in the conventional Hubbard model. In a realistic scenario for the 
Mott-Hubbard transition we expect a growing effective charge mass (or, equivalently, -Dc(O)) 
because the transport of charge becomes more difficult due to the Coulomb repulsion0. Also, 
the spin transport should become less effective because the spin exchange energy smoothly 
reduces from 0{t) to 0{J = At^ /U) and we thus also expect an increasing effective spin 
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mass (or Ds{0)). Consequently, there should be an increase of 7 below the MCIT in any 
realistic Hubbard-type model. 

The specific heat as a function of temperature is shown in fig. ^. The general structure 
of Cy{T) refiects the behavior of the density of states (see fig. |^). Even before the MCIT, 
at U = W/2, the specific heat develops a two peak structure which refiects the van-Hove 
singularity in the density of states for spin excitations. As already seen in eq. (^), the 
Sommerfeld factor is continuous at the transition. Well above the transition [U » W) 
the lower and upper Hubbard bands are well separated, which is refiected in a narrow 
low-temperature peak in c^(T) and a broad maximum around T = 0{U). In spite of the 
peculiarities of our model the overall shape of Ct,(T) is expected to be generic for any Hubbard 
model with a smooth dispersion relation, in the absence of perfect nesting. 



2. Compressibility 

Next we discuss the isothermal compressibility k(T) = dn/dfi which is related to the 
particle number fiuctuation by K{n = 1,T) = T{{AN'f)/L = - [d"^ /dfi^) f{fi,T) 
Here, AN = N — (N). A direct calculation gives 



At=C//2 



K(n = l,T) = - die [exp (-/5e,(/C)) + exp (/5e,(/C))]"' . (33) 

71 J-n/2 



Note that the compressibility probes the system slightly away from half-filling. It is thus 
seen that not only the chargeon dispersion ec(/C) but also the spinon dispersion es(/C) enters 
the expression for the compressibility. Therefore, strong charge-spin separation, in the sense 
of a completely decoupled response of charge and spin excitations to an external force, does 
not exist even at half-filling. It is only for T — ^ that the spinous do not contribute to the 
compressibility (or holons to the magnetic susceptibility), as can be seen from eqs. ( p.5b|) , 
(|15|), and ([T6D . 

The fiuctuation of the particle density at half-filling is shown in figure |. At low tempera- 
tures and below the MCIT the fiuctuations are linear in temperature, i.e., the compressibility 
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is constant and given by (|15b| ). Above the transition the charge gap opens, and there are 
only exponentially small particle number fluctuations at low temperatures. At high tem- 
peratures the fluctuations saturate at TK{n = 1,T ^ oo) = 1/2. This is the classical value 
for spin-1/2 electrons on a lattice where on average half of the sites are doubly occupied or 
empty. 



3. Magnetic susceptibility 

The magnetic susceptibility at zero external magnetic fleld is given by x{T) = 
— (9^/ dHl) fiHo, T) . It can be directly evaluated at half-fllling in terms of both spinon 
and chargeon degrees of freedom as 



x(n = l,T) = ^ rf/C[exp(/5e,(«))+exp(-/3e,(/C)r' . (34) 

TT J-7r/2 



The result is plotted in flgure |^. At low temperatures it shows Pauli behavior which is 
strongly enhanced by the interaction, especially above the MCIT (see eqs. (|15d| ), (|16f| )). 
With increasing U /W , due to the enhanced density of spin excitations at low energies, 
the susceptibility develops a strong peak at low temperatures. At high temperatures the 
susceptibility shows Curie behavior, x(n = 1,T ^ oo) = 1/(2T), the classical value for 
spin-1/2 electrons on a lattice where on average half of the sites are singly occupied. This 
behavior of x{T) is familiar from the Hubbard model with cosine dispersion at half-flUingEl, 
and is thus a generic feature of Hubbard-type models. 



B. 1/r-tJ model 

Away from half-fllling there is no compact representation of the free energy in terms of the 
spinon and chargeon dispersion and one must use the general form in eq. (pSj), with D]c = 0, 
and P]c = exp (— /3J (tt^ — /C^) /4). Furthermore, the chemical potential now depends on 
both temperature and density. We will concentrate on a typical flUing factor, n = 0.75, 
which corresponds to a critical value of the coupling Jc/W = 0.81. 
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1. Chemical potential 



The chemical potential as a function of T for fixed n = 0.75 and various values of J/Jc 
is shown in figure ^. It is seen that fi{J,T) depends smoothly on T for T/Jc < 1. Both 
above and below the MSIT fi is only weakly temperature dependent at low temperatures, 
a behavior which is especially pronounced for larger values of J/ Jc where it remains at its 
T = value, eq. ( |26a] ), for all T/J^ < 0.5. Although fi{J,T) develops stronger T-dependence 
at low temperatures in the vicinity of the MSIT, it remains continuous as a function of J 
and shows no anomalies at J = Jr. 



2. Specific heat 

The specific heat as a function of T for fixed n = 0.75 and various values of J/Jc is 
shown in figure |^. There is just one maximum which gradually becomes broader and shifts 
to higher temperatures with increasing coupling strength J. More revealing features appear 
in the temperature dependence of the Sommerfeld coefficient 7(T) = Cy{T)/T, shown in 
fig. H, which is nothing but {ds/dT), the temperature derivative of the entropy. First of 
all, it is seen that the Luttinger Liquid relation 7(T = 0) = (7r^/3)(Ds(0) + -Dc(O)) remains 
valid (compare eqs. ([T9|), ([2T|), p^), and (|3lD). For low temperatures and just above the 



MSIT, 7(T) shows a prominent peak, refiecting the large density of states for spin excitations 
just above the spin gap. As J becomes larger, this features broadens and shifts to higher 
temperatures and is completely washed out for J/Jc ~ 2. 



3. Compressibility 

Fig. P shows the isothermal compressibility K,{n = 0.75, T) for several values of J/Jc, 
(a) below and (b) above the MSIT. For J < Jc the compressibility decreases with increas- 
ing J, especially through the suppression of the peak at T ^ O.SJc- At the transition, the 
compressibility at T = doubles because of the appearance of an extra gapless holon exci- 
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tation in the spectrum (see sec. P^I B 3|) . Just below the transition this additional density of 
states for charge excitations is already present at low but finite energies and causes a sharp in- 
crease in the slope of k{T) for low temperatures such that the compressibility for J/ = 0.99 
is actually higher than for J/Jc = 0.74 in a temperatures region around T/J^ = 0.1. Above 
the transition the additional density of states for gapless charge excitations results in a new 
low temperature peak, which broadens and shifts to higher temperatures (and is eventually 
completely suppressed) as J is further increased above the MSIT. 



4- Magnetic susceptibility 

Finally, fig. |10| shows the magnetic susceptibility x('^ = 0.75, T) for the 1/r-tJ model for 
various values of J/ Jc below and above the MSIT. The susceptibility can be cast into the 
form 



T) = ^J'' dJC {exp [2p {tJC - /x)] + 4 exp [P {tJC - /x)] + 4 exp [pj (tt^ - /C^) /4] } 



-1/2 



(35) 



As expected, for J < J^^e find a finite (Pauli) susceptibility at T = 0, eq. ( p5d]) . Evaluating 
the linear term in T from eq. (|3|) one finds x(^, T 0) = 2/( Jtt^) + [8/( JTr^)^] T + 0{T'^), 
independent of n < 1. Close to the transition, however, the temperature region over which 
this expansion is valid shrinks to zero, and the susceptibility thus seems to have a negative 
temperature gradient at low temperatures close to the MSIT (J/Jc = 0.99). Just below the 
transition we have a strong increase in the density of states for charge excitations at low 
energies which already showed up in the compressibility. Consequently, the density of states 
for low-energy spin excitations is considerably reduced. This results in a low-temperature 
dip of the magnetic susceptibility close to the MSIT. 

For J > Jc the magnetic susceptibility shows activated behavior reflecting the opening 
of the spin gap, Afig (eq. (^2]) ). It is seen that the curves for J/Jc = 0.99 and J/Jc = 1-02 
qualitatively differ from each other only in a region of very low temperatures (T ^ 0. 07 Jc)- 
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IV. CONFORMAL FIELD THEORY APPROACH AND CORRELATION 

FUNCTIONS 



Since we do not know how to express the original electron operators in terms of the 
eigenstates of our effective Hamiltonian, we cannot directly calculate any correlation func- 
tions for our model. Fortunately, away from any phase transitions, our model belongs to the 
class of Luttinger Liquids^ for which the low temperature/energy behavior is dominated 
by two gapless excitations for charge and spin with linear spectrum and different veloci- 
ties Vs, Vc (charge-spin separation). It is then natural to attempt to calculate the low energy 
behavior of various correlation functions of our model by recasting our results within the 
framework of conformal field theoryiUnl and (7-ology0'&^ techniques, both of which have 
proved extremely powerful in extracting the physics of the Luttinger Liquid fixed point. In 
this section we focus on the former approach and leave for the next section the discussion 
of S'-ology. 

It is well known that T = can often be viewed as the "critical point" in two-dimensional 
classical or one-dimensional quantum field theories: correlation functions decay algebraically 
instead of exponentially for long times and/or distances, i.e., there is no intrinsic length scale. 
The assumption that, in addition to linearizing the fermionic spectrum, conformal invariance 
also holds at low energies/temperatures restricts the behavior of the lowest order 1/L- 
corrections to the ground state energy density (E^), and the energies {Ef^±) and momenta 
(P^±) of low- lying states of system of finite size, L. In turn, this information is sufficient 
to determine the long range, long time behavior of correlation functions. [An equivalent 
approach which is based on a Landau expansion around the ground state for Bethe-Ansatz 
solvable problems gives the same resultsii]. In particular, for a one component Fermi gas 
with a linear spectrum conformal invariance implies: 



Eq — Leo 



TTCf 



(36a) 



6L 




(36b) 
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P^:^~P^='^-^{h^-h-+N^-N-)+P^ . (36c) 

Here, c is the conformal charge, v is the velocity of the the right- and left-moving elemen- 
tary excitations, are their conformal dimensions, is the momentum of the sound 
excitations in the thermodynamic limit, and N"^ are integers&B. 

A particularly simple way of determining the conformal charge from the Sommerfeld 
coefficient of the specific heat is due to Affiecki. For a one-component system with right- 
and left-moving elementary excitations 

where -D(O) is the density of states for low energies. 

The information of eqs. (|36D and ( pT]) determines the large distance x, long-time t be- 
havior of correlation functions of the (primary) fields ^h±{x,t) as foUows&H 

($..(x,t)$..(0,0)) = - . . (38) 

[X — ivt) [x + ivt) 

We are interested in the long-range behavior of the spin-spin correlation function 
C^^(r, t), the density- density correlation function C^^ir^t), and the one-particle Green's 
function Go- (a;, t), 

1 ^ 

C'^r.t) = -^((n,+,,T(t) -n,+,,i(t)) (n,,t - (39a) 
^ 1=1 



1 ^ 
^ 1=1 



^\ 2 



(39b) 



GAr,t) = E(^c,,c+,,,(t)) (39c) 

where T is the time-ordering operator. The only remaining task in computing the asymp- 
totic behavior of these correlation functions is identifying the appropriate (primary) fields 
associated with the physical operators of interest. 

In what follows we use our exact results to carry out the procedure outlined above. 
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A. 1/r-Hubbard model 



1. Conformal Charge 

The obvious generalization of ( p7D to a two-component system with charge and spin 
excitations is to replace -D(O) by (-Ds(O) + Dc{0)). As could have been expected, comparing 
with eq. (|31| ) gives c = 10. 

We note that, since this identification requires finite spin and charge velocities, it breaks 
down at the MIT. This could have been inferred from the excitation spectra for charge and 
spin (eqs. (P), (^) which are no longer linear near /C = 7r/2 but behave as oc A;" with 
a = 1/2. This is refiected in the finite size corrections to the ground state energy which 



behave as J'^/ L instead of 1/L. More precisely, from eq. (|B3D , one obtains for half- filling 



A similar breakdown of conformal invariance occurs in the context of superintegrable chiral 
A'^-state Potts mo delH for AT > 3. 

It is also worth pointing out that, as shown in appendix |B|, even away from the transition, 
there are some additive corrections to the usual formula ( |36a|) in our case. More precisely, 
we find 



-^0 ~ -^^0 = — 



(41) 



6L 

(for n = 1 and U > W one has to read instead of v^). It is easy to see that the 
corrections guarantee that there are no 1/L corrections present at t/ = where Eq = Leo is 
the exact result for all L. This requirement follows from the fact that in our case the entire 
spectrum of the kinetic energy operator T can be obtained from A^ — 1 independent spin-1/2 
SU(2) algebras in each sector with fixed total momentum Q: Tq = t{Q — -k) + W Yh=i 
[N even) . Thus, in this case the algebraic structure is much simpler than the usual Virasoro 
algebra@0; this invalidates the conventional arguments based on conformal invariance and 
justifies the absence of finite size correction in the strict t/ = limit. 
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2. Conformal dimensions 



According to eqs. (|36b| ) , ( p6c| ), in order to identify the conformal dimensions we need to 
consider those lowest-lying states with the quantum numbers of the physical excitation of 
interest. We will consider two- and one-particle excitations which determine the long-range 
behavior of the spin-spin and density-density correlation functions, and one-particle Green's 
function, respectively. 

a. n < 1: Away from the transition the lowest-lying two-spinon excitation is obtained 
from the ground state by breaking a [||]-pair at /C = — vr + A/2 and JC = —it + 3A/2, and 
placing the electrons into the S = 1, = state at the same /C points. The corresponding 
excitation energy is 

Ej;^ - El^ = J(/C = -TT + 3 A/2) = ^t;f + O (l/L^) , (42) 

corresponding to h'^ = 1 (only right-moving spinous). Furthermore, 

1 9ir 
Ptt - = ^ + 2 [(-^ + ^/2) (-^ + 3A/2)] = -h-t (43) 

where we took into accountB that each bound spin-pair contributes a momentum of tt, and 
each unbound spin at /C adds /C/2. As a result, P,°? = 0. 

rig 

The charge excitation can be treated analogously and one finds (again for n <1) 

Ef^^ -E^ = {-t){ICF - 3 A/2) + t{lCF + A/2) + J{JCf - A/2) - J{1Cf + A/2) 

= ^-l\v^\ + 0{l/L^) (44) 

1 2tt 
P^- - Po" = 2 [- - 3 A/2) + {ICf + A/2)] = -, (45) 

implying = 1, P?^ = 0. 

Finally, the one-particle excitation from the A^-particle ground state is given by the 
-|- 1-particle ground state: 



ground state (A^ -|- 1 particles): o" [ti] . . . [t 
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We then find, 

- E^N,o = KlCj. + A/2) + Y! (jiJC) - J(/C + A) 

IT- _ I o A /n ^ 



27r 



IC-A/2 

E' 

K=-n+3A/2 

+ O (l/L^) 



P. 



n 



N+1,0 



Pn,o = o [(-^ + ^/2) + C^F + A/2)] = + vr(n - 1) , 



(46) 
(47) 



and correspondingly, h% = h% = 1/4, = 7r(n — 1) = kp {k% is the electronic Fermi 
momentum) . 

h. n — 1: Below the transition (U/W < 1) the same results as for n < 1 apply. 
Above the transition {U/W > 1) we have no gapless charge excitations. This implies that 
all correlation function involving charge excitations (density-density correlation function, 
one-particle Green's function) decay exponentially because their excitation energies always 
involve the charge gap A/Xc. 

On the other hand, we have additional spin excitations. Besides the excitation with two 
right- moving spinous {hf — 1, P^ = 0) we also find the corresponding excitation with two 
left-moving spinous (hj = 1, P?^ = 0). Furthermore, we may split the two spinous to form 

fig 

the excited state 

spin-exc: a['\i] . . . [\iW 



leading to 



Eh - E^ 



-K-A/2 7r-3A/2 
/C=-7r+3A/2 A:=-7r+5A/2 



27r 

T 



— + 

4 4 



+ 



o (i/l2) 



Pnf - = ^ + 2 [(-^ + ^/2) + - A/2)] = TT 



(48) 
(49) 



This gives h+ = hj = 1/4, and P^ = %. 



24 



3. Correlation functions at large times and distances 



a. n < 1 or (n = 1 and U/W < 1): According to the results in the last subsection and 
the general formulae from conformal field theory (^) we are now in position to deduce the 
long-range behavior of correlation functions: 

C^S(x,t)r^A(^^] (50a) 
C^^(a:,t)^gf \ \ (50b) 

G^{x,t) , = , (50c) 

^{x - v^t) + i/At ^{x + \vl^\t) + i/At 

where Af is a cut-off parameter (A^ = Asgnt). 

Note that the two-particle correlation functions are of the Fermi Liquid form, with 
renormalized velocities. On the other hand, the one-particle Green's function displays Lut- 
tinger Liquid behavior involving square-root singularities rather than the conventional quasi- 
particle asymptotic form, 1/ (x — vpt). Consequently, the Fourier transformed one-particle 
Green's function shows no quasiparticle peak, i.e., there is no contribution proportional to 
S{uj — vpk) in the one-particle spectral function. 

In fact, the form of Ga{uj,k) is very interesting but rather complicated. A detailed 
analysis of its properties was given recently in refs. Nevertheless, we can already 

see by dimensional analysis that there is a step-discontinuity in the momentum distribution 
^fc,(7 = (Cfco-Cfecr) — Sr 6xp(zfcr)G(r, t = 0_) at k = kp = 7i{n — 1). Our model provides 
a novel example of a system which displays a discontinuity in Uk^a in the absence of single- 
electron like quasi-particle excitations. Such systems have been termed "free Luttinger 
Liquids"lli or "Gutzwiller Liquids"0. This unusual behavior would reflect itself, for example, 
in the dependence of the Kondo-temperature on the Kondo coupling impurity embedded in 
a "Gutzwiller Liquid"i§. 

b. n = 1 and U/W > 1: The charge-charge correlation function and the one-particle 
Green's function decay exponentially, while the spin-spin correlation function has the asymp- 
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totic behavior 

C'Hx,t)^Ai^^\ +A,{^—-\ +A, , — , _ . (51) 



1 \^ / 1 \^ exp {inx) 



x-v^tj \x+\v^\tj ^(x - v^t) (x + \v^\t) 

Due to the mixture of right- and left-moving spinous, an additional structure ("2^^?- 
oscillations" ) appears in C^^{x,t), indicating strong antiferromagnetic correlations beyond 
the MCIT. 



4- Equal-time correlation functions 

Conformal field theory only allows us to calculate the long-range (large x and t) behavior 
of correlation functions. On the other hand, the large x behavior at t = 0, involves contri- 
butions from all frequencies. From eqs. ( |5(]| ) it is already clear that the low-frequency modes 
lead to a 1/x^ decay for the equal-time spin-spin or density- density correlation function. 
What can be said about the high-frequency contributions? 

a. n < 1 or n = 1 and U/W < 1: We know that there is a sharp cut-off wave vector 
for both charge- and spin-excitations. In particular, no spin-excitations are possible above 
JCf = vr(2n — 1). The highest momentum spin excitation of the ground state is 



maximum Q spin-exc: o"[tJ,] . . . [\ l]cr' 



o . . . o 

JCf 



with Q = 7T + [{JCf — A/2) + (— tt + A/2)] /2 = rnr. This corresponds to a sharp edge in the 
correlation function as a function of momentum g at q = Q (for n < 1 orn = 1 and U/W < 1 
the correlation functions do not diverge at q = Q). As in ordinary Fermi Liquids, this reflects 
itself in long-range oscillations in real space: C^^{r,t = 0) ~ Ai/r^ + A2 cos{7mr) / r"^ . The 
same argument applies to the t = density- density correlation function which also shows 
long-range 2/ci7'-oscillationsi. Finally, the single-particle momentum distribution shows step 
discontinuities at both ends of the U = Fermi "surface" . 

b. n = 1 and U/W > 1: In this regime, the entire large distance behavior is already 
contained in the large momentum-transfer modes which are also gapless. More explicitly, 
for large r we find 
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C^^r) A^-^ + A^^ (52) 

which shows strong antiferromagnetic correlations above the MCIT. The slow oscillating 
decay of C^^(r) corresponds to a logarithmic divergence of the spin-spin correlation function 
in momentum space near g = vr. 



B. 1/r-tJ model 

1. Conformal Charge 

As expected, eq. ( PID implies c = li. Note that the formulae (^) cannot be applied 
exactly at the MSIT, J = J^. At the transition we have both gapless right-moving spin- 
excitations and gapless right- and left-moving charge-excitations with finite velocities = 
JcVr/2, ff^ = — ff' = Jc'Kn/2. Nevertheless, the correct form of the 1/ L correction at J = Jc, 
Eq — Leo {J = Jc) = — [^/{QL)] T^nJci does not contain a contribution from all these three 
gapless excitations. Rather, since a finite density of spin-excitations produces an effective 
gap for the charge excitations and vice versa, only one of the two excitations, spin or charge, 
should be taken into account at J = Jc- More precisely, the factor, nnJc in the 1/ L correction 
to the energy can be written in either of two forms, imJc = vf{J J^) + \Vc{.J Jc) \ ~ 
\Vc{J = 0)1 (for the origin of the term ff'(J = 0) see the discussion for the 1/r-Hubbard 
model, eq. (p|)), or irnJc = v^{J J^) + \Vc{.J Jt)\- other words, the correct 
value of the 1/L corrections is obtained by taking the limit from either above or below the 
transition. 



2. Conformal dimensions 

a. J < Jc: We can again use all results from the 1/r-Hubbard model in the limit of 
large U /W . For two-particle excitations we then obtain h'^ = 1, = 0; = 1, = 0. 
Note that for half-filling we have additional excitations which give hf = 1/4, P°± = it. For 
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the one-particle excitations all excitations involving charges at gaped for ra = 1, while for 
n< 1 we have h% = h% = 1/4, = 7T{n - 1) = kp. 

b. J > Jc-' Since above the transition we have no gapless spin excitations, all correla- 
tion functions involving spin excitations (spin-spin correlation function, one-particle Green's 
function) decay exponentially. On the other hand, we have additional charge excitations: 
besides the right-moving holon near /Ci {h^ = 1, = 0) we also find a left-moving holon 
near /C2 {h^ = 1, = 0). Furthermore, we can have a left-moving holon near ICi and a 
right-moving holon near /C2- These excitations obey hf = 1/2, P5± = ±(/C2 — A^i) = ±2?™. 



3. Correlation functions at large times and distances 



For J < JcWe obtain the same results as for the Hubbard model away from the transition, 
while J > Jc both the spin-spin correlation function and the one-particle Green's function 
decay exponentially. On the other hand, the charge-charge correlation function behaves 
asymptotically like {vc = = —v^ = Jvm/2) 



X — Vrt 



+ 



X + Vet 



exp {i2TTnx) exp {—i27mx) 
(x - Vet) 



[x + Vet) 



(53) 



Note that since in this case large momentum-transfer excitations are also gapless, C^^{x,t) 
displays "4/cp'-oscillations"i. 



4- Equal-time correlation functions 

For J < Je we know that the exact wave function is the Gutzwiller-projected Fermi-sea, 
in which case all ground state correlation functions are explicitly knownil0 for all densities n 
and distances r. In particulaJiil, 

C^^{r ^ 0) = [Si (vrr) - Si {irr (1 - n))] (54) 

vrr 

where Si(x) = dtsint/t is the sine- integral. At half- filling, the asymptotic behavior^, 
C^^{r > 5) ~ (— l)'"/(2r), reflects the strong antiferromagnetic correlations. The formulae 
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for the density-density correlation function are rather involved but they show the expected 
large-distance behavior 



The momentum distribution is found to have a jump discontinuity's of size ^fc^^ — ''^fc|,_ = 
a/1 — n. It is comforting that all general considerations from conformal field theory are 
verified in a limit where a complete description of the ground state properties is available. 

For J > Jc, both spin-spin correlation function and the one-particle Green's function 
decay exponentially. The large momentum transfer processes for the density-density corre- 
lation function are no longer gaped and the large distance behavior even at t = can be 
deduced from conformal field theory. We find 



The momentum transform of the density- density correlation function shows a logarithmic 
divergence at g = 27m. 

V. g-OLOGY APPROACH AND CHARGE TRANSPORT PROPERTIES IN THE 



At low energies/temperatures normal electronic systems in one dimension can be de- 
scribed by a continuum field theory!^. The essential idea is to linearize the electron excitation 
spectrum near the two Fermi points, resulting in left- and right-moving fermions (momen- 
tum and/or energy transfer cutoffs A^, A, must then be introduced to regularize integrals). 
These two species of fermions interact via several scattering channels characterized by cou- 
pling constants Qi, i = 1, . . . ,4 which can also be spin-dependent ("(^-ology" Hamiltonian) . 
The two large momentum transfer processes are described by gi, which parametrizes the 
scattering process which interchanges right- and left-moving particles ( "backscattering" ) , 
and g^, which represents the scattering of two left-moving particles into two right-moving 




(55) 





(56) 



1/r-HUBBARD MODEL 
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particles and vice versa ("Umklapp-scattering"). The remaining processes, described by g2 
and Qi, involve a small momentum transfer, between a left-moving and a right-moving elec- 
tron ((72), and between electrons on the same branch ((74). The model with (72 and only is 
the Luttinger mo del§ which has been solved exactlySSS. 

In general, there is no recipe to link the g-ologj coupling constants to the parameters 
of a given lattice Hamiltonian without solving the lattice model exactlyQ. Moreover, this 
identification is only valid for interactions strengths smaller than the cutoff. In the case of 
the Hubbard model with cosine dispersion the Mott-Hubbard transition happens at half- 
filling for U = 0+0, so that the entire low-energy physics, including the physics of the 
transition, can be described within (yf-ologyEl. More generally, this approach is applicable 
to a wide class of one-dimensional electron systems, the so-called Luttinger Liquids, the 
low energy behavior of which is controlled by a weak coupling fixed pointlli. Below we will 
restrict ourselves to the discussion of the metallic phase (ra < 1, or = 1 and U/W < 1) 
of the 1/r-Hubbard model, although a similar treatment can be also given for the metallic 
phase of the t — J model. 



A. Identification of tiie g-ology parameters 

The 1 /r-Hubbard model is particularly simple as it describes only right-moving electrons 
which are, however, hopping on a lattice. It is thus immediately evident that, in the context 
of (7-ology, the corresponding low-energy physics is described by a pure g^-model ("chiral 
Luttinger model"), with a Hamiltonian 

Hg^ = Y.i^'^F^) O'ta^Ka + \Y.\Y. 9T' Pa{q)pa'{-q) ■ (57) 
fc.cr q u,u' 

Here, Po-(g) = I]fc '^k-q,a^k,ai system volume (length of the ring) is = La, Tivp = ta, and 
all terms of the Hamiltonian are understood to be normal ordered with respect to the ground 
state of the non-interacting Hamiltonian, Hq = Hg^=Q. The coupling matrix in ( pTl) can be 
decomposed as g^'" = gldcry + g't^a,-^'- For sufficiently small values of U {U/W << 1) the 
lattice plays no role, and we may obviously identify gl = (9(f/^), g^ = U . 

30 



After bosonization@'0ii the Hamiltonian (^) becomes diagonal in the new bosonic 
operators for charge (dfc) and spin {(3k) and thus, the pure (74-model is a "non-interacting" 
Luttinger Liquid or a "Gutzwiller Liquid"0. As discussed above, such a model does have a 
jump discontinuity in the momentum distribution although it is not a Fermi Liquid. The 
bosonic version of the 5f4-Hamiltonian now reads 

H',, = T.Mat<^k (58a) 

k 

Hl = T.M^Ptk. (58b) 

k 

with 

9l 



\Uc\ = 


vf 


1 + 


\Us\ = 


vf 


1 + 


9l = 


gl + gi 


9l = 


gl - 


-gi 


= U 




-gl, 



2TihvFa ^ 

gl 



27ihvFa 



(58c) 

(58d) 
(58e) 
(58f) 

I = vf{1 + U/W), \us\ = vf{1 - 
Note that the lattice provides the natural cut-off parameters, A = W and A^ = 

7r/a. 

It is amusing that, for our model, there is a way of extending the g-o\ogy solution beyond 
the perturbative regime. In particular, by identifying the velocities, |mc,s|, with the velocities 
of the charge and spin excitations found in our exact solution, eqs. (|13D, ([7a|) , we can extract 
parameters g'^^ for all values of U/W {U/W < 1 for ri = 1). The resulting expressions read 
{n = N/La): 

W 

gl = U I > (59a) 

^{W + Uf - mUna 

For na = 1 (half-filling) (74 ^ 00 for f/ W-, reflecting the MCIT in the exact solution. The 
identification (|59|), for which we give a different argument in appendix |C| (for n = 1), allows 
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us to follow the solution all the way to the MCIT by using the (/-ology parametrization. 
This clarifies the fact that, in our model, the MCIT is not caused by the the relevance of 
(73-processes at half-filling as in the Hubbard model with cosine dispersion (perfect nesting 
property), but rather, it arises as a pure renormalization of gl due to the presence of the 
lattice. 



B. Conductivity at zero temperature: Drude weight 

In this section we make use of the g-ologj approach to calculate the zero temperature 
Drude weight, Dc, of the zero-frequency peak of the real part of the conductivity. 

Re [a{uj)] = VcS{uj) for to 0. (60) 

In the non-interacting limit, Vq = e^to? / [tx?!?). 

To obtain the electrical conductivity we start with the charge-charge Green's function 
i.pi.'i) = PtI*?) + Pii'i) is density operator in momentum space), 

X^''(g,t) = ^(Tp(g,t)p(-g,0)) (61) 

the retarded part of which is easily obtained from the (7-ology analysiS@0ii: 



■K uj — v^q + irj 



This form becomes exact for our model for small values of q and tu, and implies that the 
low-energy charge transport is entirely dominated by holons, eq. (0). As a check one may 
calculate the zero temperature compressibility as k{T = 0) = limg_^o lim<^^o Xrefl?! i^) = 
1/ (vrlf^l^ in agreement with our direct calculation, eq. (|15b| ). 

Up to constant prefactors which we will put back in the end, the electrical conductivity 
follows from the generalized Einstein relation, 

aret(g,^) = ^Xref (9'^)' (63) 
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which can be derived using the hnear response formula for (JrctiQ-,^) and the continuity 
equation, —qj{q,u) = Ljh{q,u) (see also appendix 0). The real part of the conductivity can 
then be calculated from Re [aret{(l,^^)] = \Vc\S{^^ — v^q), 

cr{uj) = lim Re [(Ti.et(5', ^^)] = l'^^fl'^(i^) • (64) 

Adding back the constant prefactors gives then the final expression for the Drude weight 

U 



V, = Do— = Vo 
Vf 



1 + 



(65) 



^J{W + Uf - AWUna 
Note that in our model the Drude weight increases with interaction. In fact, at half- 
filling, Vc is seen to diverge as (Tdc(i^) = o'o(^^) / ~U /W) for small uo and U /W < 1. At the 
same time, the compressibility goes to zero, in such a way that the product, kV^, remains 
constant. This behavior is in contrast to the MCIT observed in the Hubbard model with 
nearest-neighbor hopping, where, at finite U the MCIT can be approached from the metallic 
state by increasing the filling, n ^ 1^, U > 00. In that case, the compressibility diverges 
and the conductivity tends to zero. 

In fact, it is clear that for any realistic dispersion the conductivity is bounded from above 
("f-sum rule"^iii) and cannot diverge: duKea^uj) < F < oo. As shown in appendix 0, 
F is finite whenever the bare dispersion e{k) has a finite second derivative with respect to k 
everywhere. It is thus seen that it is our discontinuous dispersion relation e{k) = t{k mod27r) 
which allows for a diverging conductivity. 

Finally, we note that since we do not know how a physical vector potential couples to 
the effective particles of our exact solution one cannot reliably determine the Drude weight 
by twisting the boundary conditioners, corresponding to adding a flux, $, through the 
ring. In fact, if we naively replace all values /C by /C -|- $/ (La) to determine the ground state 
energy Eq{^), and try to evaluate 



V = 



we would obtain the incorrect result for the Drude weight 

33 



(66) 

$=0 



V = Vo 



— - 1 + 

Vp Vp 



(67) 



which involves contributions from both charge and spin excitations. 



VI. SUMMARY AND CONCLUSION 

Above we discussed various aspects of the physics of two exactly solvable models of one 
dimensional lattice fermions: the 1/r-Hubbard model and a related 1/r-tJ model. The 
special feature of these models is that, due to absence of perfect nesting, they display non- 
trivial metal-to-charge-insulating and spin-insulating states, respectively, at a finite value of 
the coupling constants. 

The 1/r-Hubbard model thus provides for an explicit realization of Mott's and Hubbard's 
ideasi'i albeit in a rather pathological model in which the dispersion has a jump discontinuity 
at the zone boundaries. This peculiarity allows for a diverging Drude weight when the 
MCIT is approach from the metallic side, associated with the vanishing of the charge mass 
at the transition. While the behavior in the charge sector in the 1/r-Hubbard model is 
rather non-generic, the spin sector reflects the expected physics of correlated Fermi systems: 
for example, the magnetic susceptibility strongly increases with interaction and develops a 
strong low-temperature peak at the spin-exchange energy scale J = At'^/U. After the MCIT 
the static spin-spin correlation function displays strong antiferromagnetic correlations with 
a logarithmic divergence in momentum space at g = vr . 

The 1/r-tJ model with pair-hopping terms displays a metal-to-spin-insulator transi- 
tion at Jc = 4t/(7r(l — n)). The spin liquid state is described by the Gutzwiller projected 
Fermi-sea and provides an exactly solvable example of a genuine RVB-state. At the MSIT, si- 
multaneously with the opening of a spin gap, the zero-frequency "spin-conductivity" (Drude 
weight for spin transport) drops to zero, after having increased with J/t (to a finite value) 
on the spin liquid side. The MSIT appears as a pure level-crossing effect, and thus should 
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not be interpreted as a Mott-type transitionQ. 

We also found that, with a few changes related to the high symmetry of the spectrum in 
the non-interacting {U = and J = 0) models, conformal field theory is applicable in the 
metallic and spin liquid states (i.e., away from transitions). Using conformal field theory 
and g-o\ogy techniques allowed us to calculate long-distance/time properties of ground state 
correlation functions. The two-particle sector shows Fermi Liquid behavior. However, the 
single-particle Green's function does not lead to a quasi-particle peak in the spectral function, 
even though the jump discontinuity in the momentum distribution survives. Such systems 
are referred to as "free Luttinger Liquids" or "Gutzwiller Liquids". Within the g-o\ogy 
approach the 1/r-Hubbard model is identifyed as a "pure g^-model" or "chiral Luttinger 
model" . In our model the MCIT is driven by short distance lattice effects which lead to a 
divergence of the interaction parameter gl at the MCIT. 

Despite their conceptual shortcomings the two models provide instructive examples for 
strongly correlated electron systems. Although a rigorous proof is still missing we think that 
our models and results can be used to check the capability of numerical techniques and the 
applicability of various approximation^!. 
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APPENDIX A: CALCULATION OF GROUND STATE PROPERTIES IN THE 

PRESENCE OF A MAGNETIC FIELD 



In this appendix we calculate some ground state properties in the presence of a weak 
magnetic field Ho, and will eventually let Ho 0. 



1. 1/r-Hubbard Model 

In the presence of a magnetic field we will have a finite magnetization. Spin pairs will be 
broken and turned into the direction of the magnetic field. The positions of these unpaired 
spins is determined by the spinon dispersion (||). From the discussion above and figure |I| it 
is clear that they will always be located around /C = — vr. For n = 1 and U/W > 1 there are 
also some upturned spins around /C = vr. We exemplify the calculation for the latter case. 

For n = 1, U/W > 1, and m > the ground state is represented by 



ground state: T • • • T 



[Ti]---[n] 



T...T 



ICo 



where the magnetization m is given by /Ci — /C2 = 27r(m — 1). The ground state energy 
density is calculated as eo{n = 1,/Ci,/C2) = {l/2n) Jl^^ die [J {IC)/2] where the additional 
factor 1/2 took into account that only every second /C value contributes to the integral. We 
are only interested in small fields and expand /Ci,2 = =F7r(l — ?7i_2) with r^i ^ ^ 1. To second 
order we then obtain with the help of eqs. @, (^), and (|7bD 

eo{vi,V2,n = l) = eoin = l) + ^{r]fvf + r]j\v^\) . (Al) 

Using m = {rji + f]2)/2 we minimize 60(71 = 1, 771, 772) — mHo with respect to r^i 2, and obtain 
Vi,2 = (27^o)/(7r|'yf'^|)- The magnetization becomes 

™ = H„(K))-' + (.K'|)-')=l^ (A2) 

and the ground state energy reads 

1 W'^ 

eo{m ^ 0, n = 1) = eo(n = 1) + --rrrn? . (A3) 

o U 
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A similar calculation gives the ground state energy density for n < 1 or n = 1, U/W < 1 



as 



eo{m, n < 1) 



U{n — m) — W [{1 — n)n + (1 — m)m] 



2AWU 



({W + Uf - mUmf^ - ({W + Uf - AWUn 



3/2 



(A4) 



which is valid for all m. Note that the limits TYq and n ^ 1 do not commute beyond 
the metal-to-insulator transition, just as in the usual Hubbard model@. 



2. 1/r-tJ Model 

The ground state representation in the presence of a (not too big) magnetic field is given 



by 



ground state: 





• T 


o . 


. o 











/Ci 



[Ti] • • • [Ti] 



1C2 



with /C2 — /Ci = 271(71 — m) and Km = vr(2m — 1). A straightforward integration gives 
eo(/Ci, IC2, ICm) = (-t/(47r)) [JCj - JC^ + vr^ - /Cf] - ( J/(167r)) [n^iJC^ - /Ci) - (/C^ - /C?)/3]. 
We have to minimize eo(/Ci, IC2, /Cm) — fin — Hom with respect to /Ci, IC2, /Cm, where we have 
to distinguish the cases (i) /Ci = /Cm, and (ii) /Ci > /Cm to allow for a transition at some 

J=Jc- 

(i) /Ci = ICm- we are below the transition, and we see that /C2 = 7r(2r2 — 1), /Cm = 
7r(2m — 1). In this case we can find the ground state energy density even without the 
minimization procedure. For small external fields it has the simple form 

Wn{l — n) J7r^(n — m) 



eo(n,m, J < Jc) 



3(n + m)— 2(n + nm + m ) . (A5) 



2 12 

By differentiating this expression with respect to m gives the connection between the external 
magnetic field ?io and the (small) magnetization density m 



m = Ho (nv 



.R 



-1 



(A6) 
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(ii) /Ci > ICm'- from the minimization equations one easily finds that /Ci + /C2 = 
—AW/lnJ). We thus get a solution 

2W 

/Ci > /C^ , if J > J,{m) = —- . (A7) 

TT^[l — n — m) 

Note that the presence of a magnetic field stabilizes the spin-liquid phase where spin- 
excitations are gapless: Jc{m) > Jc{m = 0) = Jc- The ground state energy density then 
reads 

eo{n,m, J > Jc{m)) = — ~ '"^) m(l — m) {n — ni) 3 — {n — m) . (A8) 

By differentiating this expression with respect to m gives the connection between the external 
magnetic field Tio and the (small) magnetization density m 

no = n^o + rmr^ [J^ + n{J - J,)] (A9a) 
n^o = ^^^^^^{^-j)W + n)-Jc{l-n)] . (A9b) 

It now takes a finite magnetic field Tip to magnetize the system. The spin gap is A/i^ = 27ig 
(see (13)) because we break a pair of spins {S = 0) to form an = S = 1 state. 



APPENDIX B: 1/L CORRECTIONS FOR THE GROUND STATE ENERGY OF 

THE 1/r-HUBBARD MODEL 

The ground state energy for finite system sizes L (even) and a finite particle number N 
(even) is calculated for < 1 as Eq = -tt<k<Kp ^ic + J2kf<ic<tt{^'I^)^- The prime 
indicates that only every second of the /C = A(m + 1/2) m = —L/2, ... , L/2 — 1) has to 
be taken {JCp = A{N — L/2)). After rescaling m = 2r — L/2 + 1 the sums give 

Un Wn{l — n 



E- = L 



^ N/2-1 

^ E \/W'^ + U^ -4WU{2r+l)/L . (Bl) 



4 4 

To calculate the sum we use the Poisson sum formula (see, e.g., Ref. ^9|) 



V/i(r + l/2)=/ c/x/i(x) + 2 V(-l)' / dxh{x) cos {27TXS) , (B2) 

r=0 .=1 
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do the first integral and a partial integration in the second, and arrive at 

r ^ ,^s2W [N/^ sin(27ra;s) 

Ltis Jo ^I{W + Uf -mUx/L 



(B3) 



Away from U = W and n = 1 we can do further and further partial integrations in the 
integral on the right hand side of this equation, and generate an expansion in 1/L. Keeping 
only the first term gives 



E^-Lto = -Y. 



WU 



\ s2 L7r2(iy + f/) 



1 



1 



. (B4) 



'1 -4Wn/(W^ + f/) . 
Doing the remaining sum and rearranging terms {vp = t = W/{27t) is the Fermi velocity) 
gives 

U U 



Eq - Leo = TrrVp 
bL 



U + W ^(^w + Uy ~AWUn_ 
With the help of eqs. (0) and (|TD|) this can be cast into the form of eq. (PT|). 



(B5) 



APPENDIX C: EFFECTIVE g4-PARAMETERS AT HALF-FILLING FROM THE 
SCREENED ELECTRON-ELECTRON INTERACTION 

There is a rather simple recipe to obtain gl'^ for all U/W < 1 without referring to the 
exact solution. At half-filling we know that kp = lies symmetrically around the upper and 
lower band edge {±W/2 aX k = ±71/ a). For a pure (74- model one can exactly calculate the 
screened electron-electron interaction with the help of Ward-identitieJ^^Sil. One finds 

u — Vpk 



Dl''{k,uj) = gl'' 



(CI) 



u — \Uc,s\k 

The effective interaction integrated over all times is given by the 00 = contribution which 
gives 

Vp . » 1 



Dr{k,o) = gr 



c,s 
94 



(C2) 



i+gr/w • 

If we now demand that the time-integrated part of the screened interaction has to stay 
unrenormalized at its value for small U/W, Dl{k, 0) = U and Dl{k, 0) = —U, we find that 
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which is indeed the correct result. So far, we have not been able to obtain a general form 
for 514* away from half-filling, without appealing to the exact solution. 

APPENDIX D: f-SUM RULE FOR GENERAL DISPERSION RELATIONS 

1. Current operator 

For a Hamiltonian H — T -\-V v^e assume that the interaction only contains density 
operators fir. We omit spin indices in the following. The particle density operator n(g) = 
^k+q^k then commutes with V . We expand the exponential in the particle 
density operator Ti{q) for small momenta q (this step could cause problems in the case of 
long-range hopping). The continuity equation then reads 

Q 

i-Q^n{q) = -q'^—fir^ -q'^v{r)n{r) . (Dl) 

The sum on the right-hand side is identified with the particle current operator 

^ — % d 

7 = lim — 7— n((/) . (D2) 
</-o q dt ^ ' ^ ' 

Using the Heisenberg equations of motion we can calculate j as 

If the dispersion relation is differentiable, this reduces to the standard expression for the 
current operator involving the group velocity de{k)/dk. Our linear dispersion relation is 
not differentiable at \k\ = tt. Instead we may think that the linear dispersion relation is 
the result of a limiting process where we start with a differentiable dispersion e{k) which 
becomes ever steeper near — tt. This corresponds to ever faster left- moving electrons 
near = tt, in addition to the right-moving electrons for |A;| < tt. From this viewpoint 
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it becomes clear that we have to expect a singular contribution from these "left-moving" 
electrons. Indeed, we find for g > 



W 



2^ ^ ^k+q(^k ~ ^k+q^k 



w ^ 



(D4) 



^ -^^-7r+A/2^7r-A/2 



and j{—q) = Now let g — and assume that we may approximately replace c^^_q 

c^. For the minimum q = A = 2tx / L one has j{q = A) = W/{2tx) 
As "the current at g = 0" one may define 
^ 1 



j = -(j(g = A)+j(g = -A) 



l^/(2vr) [n - (L/2) (cl^^^/2C._A/2 + ^:_a/2^-.+a/2) 



(D5) 



We obviously get an extensive (i.e., singular) contribution from the states near the Brillouin 
zone boundary. 

2. f-sum rule 



We are interested in the real part of the transverse conductivity for small q at temperature 
T = 

^2 



Re \(j{uj, g — 0)] = - — lim Re 

Luj q^o 



dt'e^-^'-''\j^{q,t),j{q,t') ) 



(D6) 



Using the continuity equation ( |D1|) , the equation of motion (in t) (eq. (P3|)), and integrating 
by parts (in t') one finds 

^2 



Re [ct(u;, g 0)] 



Lujq^ 



Re (-z)( /3+(g,t),f ,p(g,t') )e-(*-*') 



(D7) 



The first term in {. . .} is purely imaginary and drops out. To obtain the f-sum rule we 
integrate over all uj. This gives 



F(g) 



Tie 



P+(g),Tl,p(g)l) 



(D8) 
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because the expression ( P7|) only depended on {t — t'). The double commutator gives 



P'^{(1),t] ,p{q)] = J2^tck 

k 



e{k) - e{k + g) + e{k) - e{k - q) 



Thus, 



2 ri 

k k1 



e{k + q)-e{k) t{k - q) ~ t{k) 



{fik) 



q q 

and F = limg^o -^(o')- For twice differentiable dispersion relations this reduces to F 



(D9) 



' / L)J2k{f''k) [d'^e{k)/dk'^]. In particular, for a cosine dispersion, F = 7re^(— T)/L which 



is the standard result^S. 

For the linear dispersion we get 

Tre^W 



F{q) 



Lq^ 



— 7r</i,'<— 7r+g 

The minimum value of g is A = 27r/L. Thus 

^2 



7v—q<.k<7T 



(DIO) 



F = F{A) = L— [(n,_A/2) - {n-n+A/2) 



(Dll) 



The f-sum rule is thus extensive and leaves room for a diverging Drude weight. This can 
only happen, if the dispersion relation e{k) is not twice differentiable with respect to k. 
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FIGURES 

FIG. 1. Dispersion relations for (a) spinons, see eq. @, and (b) chargeons, see eq. (P), in the 
1/r-Hubbard model for (a) U/W = 0.25, 0.5, 0.75, 1, 2, and (b) U/W = 0.5, 1, 2. 

FIG. 2. Total density of states (a) Dg of the spinons and (b) Dc of the chargeons for the 
values of fig. || in the 1/r-Hubbard model (note the change in energy scale); the energy gap 
in Dc for U > W is only half as big as the true gap because the chargeons always come in pairs. 

FIG. 3. Specific heat as a function of temperature for the half-filled 1/r-Hubbard model 
for various values of U/W (a) below {U/W = 0,0.5,0.75,1), (b) above (U/W = 1,2,4) the 
metal-to-charge- insulator transition. 

FIG. 4. Fluctuation of the particle density kBTK{n = 1,T) for the half-filled 1/r-Hubbard 
model as a function of temperature for U/W = 0, 0.5, 1, 2. 

FIG. 5. Magnetic susceptibility = 1; T) for the half-filled 1/r-Hubbard model as a function 
of temperature for U /W = 0, 0.5, 1, 2. 

FIG. 6. Chemical potential /u(n = 0.75, T) in the 1/r-tJ model as a function of temperature 
for J/Jc = 0.5,0.99, 1.02, 1.2, 2; the metal-to-spin-insulator transition happens at Jc = 0.81VK. 

FIG. 7. Specific heat c^,(n = 0.75, T) in the 1/r-tJ model as a function of temperature for 
J/Jc = 0.49,0.99,1.02,2.47. 

FIG. 8. Sommerfeld coefficient 7(77. = 0.75, T) in the 1/r-tJ model as a function of temperature 
for J/Jc = 0.49, 0.99, 1.02, 2.47. 

FIG. 9. Compressibility K{n = 0.75, T) in the 1/r-tJ model as a function of temperature (a) 
below (J/Jc = 0.49,0.74,0.99) and (b) above (J/Jc = 1.02, 1.23, 1.73,2.47) the MSIT. 

FIG. 10. Magnetic susceptibility xi^ = 0.75, T) in the 1/r-tJ model as a function of tempera- 
ture for J/Jc = 0.49, 0.74, 0.99, 1.02, 1.23, 2.47. 
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